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Target  Patterns  in  Reaction- Diffusion  Systems 

Patrick  S.  Hagan*  + 

Department  of  Mathematics.  Stanford  University.  Stanford,  California  94305 


We  consider  the  general  reaction-diffusion  system  A,  =  F(A)  +  t  £>w  v  2A + 
<g(x.  A),  0  <  (  «  I,  where  the  small  term  <g( x.  A)  represents  the  effects  of  localized 
impurities.  We  assume  that  the  system  A,  =  F(A)  has  a  stable  time-periodic  solution. 
Then  we  construct  stable  target  pattern  solutions  of  the  full  system.  For  typical 
initial  conditions  we  find  that  these  target  patterns  will  arise  only  if  g(x.  A)  ^  0. 
Finally,  we  determine  how  target  patterns  interact  and  show  that  higher  frequency 
target  patterns  eventually  engulf  neighboring  lower  frequency  target  patterns. 


1.  Introduction 

In  some  chemical  and  biological  systems  confined  to  two  spatial  dimen¬ 
sions,  target  patterns  are  commonly  observed:  concentric  circular  concentra¬ 
tion  waves  which  expand  outward  with  new  waves  being  generated  at  the 
center.  For  example,  Belousov-Zhabotinski  (B-Z)  reactions  [1-4]  exhibit 
target  patterns,  as  do  fields  of  aggregating  D.  discoideum  [5-8].  (See  Fig.  1.) 

Target  patterns  have  been  investigated  previously  [9-15].  Here  we  con¬ 
sider  the  general  system 

A,  =  F(A)  +(DmV2A  +  eg(£,A),  0  <  e  «  1 .  (1.1) 

In  (1.1)  the  diffusion  term  is  small  because  we  have  chosen  to  look  for 
solutions  which  vary  on  long  spatial  scales.  We  have  done  this  since 
diffusion  tends  to  make  A  spatially  uniform.  Also,  (1.1)  is  spatially  inhomo¬ 
geneous  due  to  the  small  term  cg(f,  A).  For  the  B-Z  reaction,  this  term  can 
arise  from  particle  contaminants  and  other  impurities  [3].  For  fields  of 
D.  discoideum ,  this  term  arises  from  the  different  starvation  states  of  the 
individual  amoebas  [6]. 

We  will  construct  formally  stable  target  pattern  solutions  of  (1.1),  first  for 
the  case  g(x.  A)  =  0  and  then  for  the  general  case  g(x,  A)  ^  0.  We  will  find 
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J^IG.  I  (Photograph  reproduced  with  permission  from  reference  (5]).  Photograph  of  the  amoebas  DUtyoswlium  Jiscoideum  in  the  aggregation  phase, 
light  bands  are  bands  of  moving  cells.  Both  target  patterns  and  spiral  waves  are  present. 
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that  the  target  patterns  for  g(x,  A)  =  0  occur  only  for  special  initial  condi¬ 
tions.  However,  for  g(x,  A)  ^  0  they  arise  even  if  the  system  is  initially  well 
stirred. 

This  matches  the  experimental  observations  of  the  B-Z  reaction  in  [3J. 
There  it  is  reported  that  there  are  small  particles  at  the  centers  of  most 
target  patterns,  and  that  filtration  of  the  chemicals  to  remove  impurities 
eliminates  nearly  all  of  the  target  patterns. 

We  also  determine  how  target  patterns  interact.  We  find  that  higher 
frequency  target  patterns  engulf  nearby  lower  frequency  patterns.  This 
matches  the  observations  in  [1]. 


2.  The  Governing  Equation 

We  assume  that  at  c  =  0,  system  (1.1)  has  a  stable  time-periodic  solution 
A(r,  x)  =  B(r)  =  B(f  +  P)  with  period  P.  We  now  use  two-timing  to  solve 
(1.1).  We  let  T  =  a,  we  expand  A  as 

A(«,  t,  x )  =  A°(r,  T,  x )  +  eA'(t,  T,  x)  +  •  •  • ,  (2.1) 

and  we  require  A1,  A2, ...,  to  be  bounded  in  t.  Substituting  (2.1)  into  (1.1) 
we  find 

A0,  -  F(A°)  =  0,  (2.2) 

A',  -  Fa (A0 )A‘  =  -  A°r  +  Dm  v  2A°  +  g(  x.  A0 ) ,  (2.3) 

Here  the  matrix  /^(A0)  is  the  derivative  of  F(A)  with  respect  to  A. 

The  solution  of  (2.2)  of  interest  is 

A°  =  B(f +  «Kr,x)),  (2-4) 

where  the  function  \p(T,x)  is  undetermined  at  this  stage.  Substituting  (2.4) 
into  (2.3),  we  obtain 

A1,  -  ^(B)A1  =  -B'*r+  DmB'V2*  +  Dm B"V*  •  V*  +  g(x,B), 

(2.5) 

where  the  argument  of  B,  B',  and  B"  is  /  +  \p(T,  x). 

Now  the  equation 


u,  -  FA(B{t  +  *))u  =  0 


(2.6) 
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has  a  periodic  solution  u  =  B\t  +  ^).  To  ensure  that  A0  =  B(/  +  is  a 
stable  solution  of  (2.2),  we  assume  that  all  solutions  of  (2.6)  which  are 
linearly  independent  of  B'(*  +  $)  decay  exponentially.  Under  these  condi¬ 
tions  there  is  a  unique  row  vector  zT(t  +  \p)  which  satisfies 

zf  +  zrF/<(B[/  +  ^])  =0r,  zr(t  +  i)B'(t  +  xf)  =  1,  for  all/, 

(2.7) 

and  is  periodic  with  period  P.  Moreover,  the  solution  A '(/.  T,  x )  of  (2.5)  is 
bounded  in  /  if  and  only  if 

fPzT[-B’4,T+  DMB'v2*  +  DMB"X>t-  Vi|/  +  g(x, B)|  ds  =  0,  (2.8) 
■'o 

where  the  argument  of  zr  and  B  is  s  +  ^(T,  x).  Using  periodicity,  we 
rewrite  (2.8)  as 

D(v2+ +  Tvt  ■  V$)  +  a(x),  (2.9) 

where  ;he  constants  D,  T  and  the  function  a(x)  are 

D  =  fPzT(s)DMB’(s)ds/P ,  DT  =  f'>zT(s)DMB"(s)  ds/P , 

Jo  Jo 

a(x)  =  fPzT(s)g(x,B(s))  ds/P.  (2.10) 

■'o 

In  summary, 

A(/,  x)  =  B(/  +  ip(T,  a))  +  0(e),  (2.11) 

where  \p  evolves  on  the  long  time  scale  T  =  ef  according  to  (2.9).  Equation 
(2.9)  is  the  key  equation  that  describes  how  target  patterns  evolve  and 
interact.  The  derivation  of  (2.9)  from  (1.1)  is  very  similar  to  a  result  of  John 
Neu  (16J. 

To  interpret  the  a(x)  term  in  (2.9),  consider  (1.1)  with  the  diffusion 
matrix  set  to  0.  Since  e  «  1,  at  each  spatial  position  x  there  must  be  a 
stable  time-periodic  solution  near  B (/).  To  find  this  solution,  we  note  that 
Dm  =  Q  in  (2.10)  implies  that  D  =  0,  so  (2.9)  shows  that  t]/  =  a(x)T. 
Substituting  this  into  (2.11),  we  find  that  A  =  B(t[l  +  ta(x)  +  ••'])  + 
0(e).  Thus  the  a(x)  term  in  (2.9)  represents  the  local  first-order  frequency 
shift  due  to  the  inhomogeneous  term  eg(x,  A). 
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3.  Spatially  Homogeneous  Systems 

Here  we  consider  systems  with  g(x,A)  =  0.  Then  (2.10)  shows  that 
a(x)  =  0,  so  (2.9)  becomes 

+T  =  d[v2^  +  rvij-  •  v>^].  (3.i) 

To  solve  (3.1)  we  apply  the  Cole-Hopf  transformation 

Z  =  er+.  (3.2) 

finding 

ZT=DV2Z.  (3.3) 

The  initial  value  problem  for  (3.1)  can  be  solved  by  solving  the  initial 
value  problem  for  (3.3).  Thus,  in  n  space  dimensions 

ip(7\.x)  =  r  '  log  |(4ttD7'  )  "/2Jexp[r^(0,  v) 

-  ( x  - y )  •  ( x  ~y) / ADT ]  dy J  .  (3.4) 

Now  (3.4)  shows  that  if  V'fO,  x)  is  bounded,  then  \p(T,  x)  goes  to  a  constant 
as  T  becomes  large.  Thus  when  g(x.  A)  =  0,  no  target  patterns  evolve  for 
bounded  initial  conditions. 

Next,  we  construct  stable  target  pattern  solutions  of  (1.1)  with  g(x.  A)  =  0. 
These  solutions  must  have  ^(0,  x)  unbounded,  and  so  they  may  be  of 
mathematical  interest  only. 

By  solving  (3.3)  and  using  (3.2),  we  find  that  (3.1)  has  the  solutions 

*(7,x)  -W7/T  +  T  %g/0(vVZ)r),  (3.5) 

in  two  dimensions  and 

i /'(T,  x)  =  uT/T  +  F  " 1  log{r^  1  sinhyto/Dr) ,  (3.6) 

in  three.  Here  r  =|  x\ ,  /0  is  the  zero-order  modified  Bessel  function,  and  w 
is  any  positive  constant.  Moreover  (3.4)  shows  that  the  solutions  in  (3.5)  and 
(3.6)  are  stable  to  all  bounded  initial  perturbations. 

Equations  (2.11)  and  (3.5)  yield  target  pattern  solutions  in  two  dimen¬ 
sions.  Similarly,  (2.11)  and  (3.6)  yield  time-periodic,  spherically  symmetric 
solutions  in  three  dimensions.  When  r  »  1,  both  (2.11),  (3.5)  and  (2.11), 
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(3.6)  reduce  to 

A(/,  x )  ~  B([l  +  tu}/T]t  +  i^/Dr/Y  +•••)+  0(c).  (3.7) 

Since  c  «  1,  these  waves  are  ingoing  if  T  >  0  and  outgoing  if  T  <  0. 

Target  patterns  are  observed  in  the  B-Z  reaction  and  in  D.  discoideum, 
even  if  the  systems  are  well  stirred  initially  so  that  4/(0,  x)  is  nearly  constant. 
In  the  next  section  we  show  that  the  observed  target  patterns  can  be 
explained  if  the  inhomogeneous  term  g(x,  A)  is  not  identically  zero. 


4.  Spatially  Inhomogeneous  Systems 

Here  we  consider  the  inhomogeneous  case,  where  g(x,  A)  0  in  (1.1)  and 
hence  «(x)  0  in  (2.9).  Specifically,  we  assume  that  a(x )  is  smooth  and 

that  fYa(x)dx  >  0,  where  the  integration  is  over  all  space.  For  mathemati¬ 
cal  clarity,  we  also  temporarily  assume  that  a(x)  =  0  for  all  |x  —  x0|>  /?. 
for  some  fixed  x0  and  R.  With  these  restrictions  on  a(x),  we  will  now  solve 
(2.9)  for  bounded  initial  conditions  i//(0,  x).  We  will  find  that  the  solution 
4/(F,  x),  together  with  (2.1 1),  will  yield  a  single  target  pattern  centered  near 
x  =  x0.  This  target  pattern  will  depend  only  on  a(x),  D,  and  I\  and  will  be 
independent  of  the  initial  condition  \p(0 ,  x ). 

First,  we  apply  the  Cole-Hopf  transformation  (3.2)  to  (2.9),  finding 

ZT=  Dv2Z+  r«(x)Z,  Z(0.x)  =  er+(0-7>.  (4.1) 

Next,  we  separate  variables  Z(T,  x)  =  eaT$(x),  obtaining  the  eigenvalue 
problem 


Z)v24>  +  Ta(x)$  =  4>(x )  bounded.  (4.2) 

With  the  restrictions  placed  on  a(x),  it  is  known  [17,  18]  that  in  two 
dimensions  there  are  a  finite  number  m  >  I  of  discrete  real  eigenfunctions 
<t,(x),  $2(x ), . . . ,  $m(x).  Their  eigenvalues  . . . ,  um  satisfy 

w,  >  «2  >  «3  >  ■  •  •  >  wm  >  0.  (4.3) 

Also  the  eigenfunction  4>,(x)  with  the  largest  eigenvalue  w,  is  of  one  sign 
for  all  x.  Furthermore,  for  each  j  —  1 . «, 


<*>,(x)~c/0)(*/r,/2e-*' 


when  kjr  »  1, 


(4.4) 
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for  some  Cy(0).  Here 

kj=fc/D,  (4.5) 

and  r  and  6  are  defined  by  x  —  x0  =  r( cos  0,  sin  0). 

Additionally  there  is  a  two  dimensional  continuum  of  eigenfunctions 
<t >(k,  3c).  These  functions  have  the  eigenvalues  -D\k\2,  and  are  the  solu¬ 
tions  of 


x)  =  e‘^+  ^  /  H<"(\  k  1 1  x  -y\)a(y)*(k,  ?)  (4-6) 

where  H^u  is  the  zero-order  outgoing  Hankel  function.  Finally,  the  continu¬ 
ous  eigenfunctions  plus  the  m  discrete  ones  form  a  complete  set  over  the 
space  of  C2  functions. 

We  now  solve  (4.1).  First  we  define  Z(T,x)  =  H(T,x)  +  u(T,  x),  where 
H(T,x)  is  the  solution  of 


HT=Dv2H,  H(0,  x)  =  Z(0,3t)  =  (4.7) 


Then  u  is  the  solution  of 

uT=  DV2u  +  ra(x)u  +  ra(x)H(T,x),  u(  0,3t)=0.  (4.8) 

At  each  T,  a(x)H(T,  3c)  is  in  £2,  so  we  can  solve  (4.8)  by  expanding  u  and 
a(x)H(T ,  3c)  in  terms  of  the  eigenfunctions  $,(3?)  and  4>(k,  x).  This  yields 

f(T,x)  =  r~' log  Z(T,x),  (4.9) 

where 

m 

Z(T,x)  = /3,eu<T<P,(x)  +  2  +  H(T,x) 

7  =  2 

m 

-  2  Y,(F)<I>y(x)  +  v(T,  x).  (4.10) 

7=1 

Here, 

fy=f $,(*  )er*(0- x)  dx  ■  |  j $f(x )  dx 


Yj(T)  =/ <bj(x)H(T,  x)  dx  •  f<bj{x)dx 


- 1 


(4.11) 
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and  t>  is  the  contribution  from  the  continuous  spectrum 
v(T,x)  =  j t](T ,  k)$(k,  x )  dfc, 

t}(T,  k)  =  F fTe~Dkls^f  $*(k,  x)a(x)H(T  -  s,  x)  dx  ds,  (4.12) 
where  *  denotes  complex  conjugate. 

Now  \p(0, 3c)  is  bounded,  so  there  are  constants  C,  and  C2  such  that 
0  <  C,  <  H(T,x)  <  C2  for  all  x  and  all  T  >  0.  Thus  (4.1 1)  shows  that  each 
y j(T)  is  also  bounded.  In  Appendix  I  we  show  that  for  some  constant  C3, 

|  v(T,  x)  | <  C3T  for  all  x  and  all  T  >  0.  There  we  also  show  that  |  v(T,  3c )  | 
<  /J1e"|T$,(3c)  +  H(T,  x)  at  all  x  when  T»  1.  Finally,  (4.3)— (4.5)  show 
that  each  of  the  first  m  terms  in  (4.10)  decays  exponentially  in  space,  but 
grows  exponentially  in  time,  with  the  first  term  growing  at  the  fastest  rate. 
Therefore,  in  the  region  where 

eutT-kt |*- ?„!«!,  T »  1,  (4.13) 

the  H(T ,  x)  term  is  much  larger  than  all  the  other  terms  in  (4.10).  Using 
(4.9)  and  (2.1 1),  we  thus  find  that 

A(f,  x)  *  B(r  +  r-'  log  H(T,x)).  (4.14) 

Moreover,  (4.7)  shows  that  H(T,  x)  goes  to  a  constant  as  T  becomes  large. 
So  we  conclude  that  in  the  region  where  (4.13)  holds,  the  system  evolves  into 
a  spatially  constant  bulk  oscillation. 

On  the  other  hand,  in  the  region  where 

e<-,T-*,|i*-r0i3>  (4.15) 

the  first  term  in  (4.10)  is  much  larger  than  all  the  other  terms.  Noting  that 
j8|4>,(x)  >  0  for  all  x,  we  rewrite  (4.9),  (4.10)  as 

*(7\x)  =  w,r/r  +  r-'iog^.tx)  +  r-'E(7\  *).  (4.16) 

Here  the  error  term  E{T,  x)  is 
£(7\  x)  =  logjl  +  2  e~(u'~u')T%{x) 

H(T,x)-  2  Yy(T)^(je)  +  «(T,x)J  /[j8,0,(x)] 

(4.17) 
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Equations  (4.3)-(4.5)  show  that  E(T,  x)  «  1  when  (4.15)  is  satisfied.  They 
also  show  that  E(T ,  x )  decays  exponentially  in  T  for  fixed  x.  Together, 
(4.16)  and  (2.1 1)  yield  our  target  pattern  solution 

A (t,  x)  =  B([l  +  €W|/r]t  +  r_l  log  <&,(*)  +  r-’logh.)  +  0(e), 

(4.18) 


where  we  have  omitted  the  small  error  term  E(T,x). 

We  now  examine  the  qualitative  behavior  of  the  target  pattern  solution 

(4.18).  First,  the  bulk  oscillation  B(t)  is  periodic  with  period  P,  so  (4.18)  is 
time-periodic  with  period  Pt—  P/(  1  +  ewj/T).  Moreover  (4.4)  shows  that 
when  /c ,  |  x  —  x0 1  =  k{r  »  1,  (4.18)  reduces  to 

A(f,  x)  =  B([l  +  ew,/T]/  -  kxr/T  +  T~ 1  log b]Cl(0)(kir)  ,/2) 

(4.19) 


So  when  &,r  »  1,  the  target  pattern  is  asymptotic  to  a  cylindrical  wave  with 
a  radial  phase  speed 


Vph  =  r(l  +  «w,/r)/Ar,  =  sgn(r)fD\T\ 


P  P 

P.P-P, 


1/2 

.  (4.20) 


In  particular,  the  target  pattern  is  outgoing  if  T  >  0  and  ingoing  if  T  <  0, 
the  opposite  of  the  target  patterns  in  Section  3. 

Second,  the  target  pattern  solution  (4.18)  is  valid  only  in  the  region  where 
(4.15)  is  satisfied.  Since  T  =  ef,  the  region  occupied  by  the  target  pattern 
spreads  outward  from  x  =  x0  at  the  constant  radial  speed 


sp 


=  ]jtD  |  r  | 


(4.21) 


Recall  that  the  diffusion  matrix  in  (1.1)  is  c DM.  We  can  make  all 
dependence  on  t  explicit  by  replacing  x  by  'Jlxd.  In  terms  of  the  new 
variable  xrf,  the  diffusion  matrix  in  (1.1)  is  just  DM  and  the  velocities  vph  and 
vsp  in  (4.20)  and  (4.21)  are  multiplied  by  c~,/2. 

Finally  except  for  the  overall  phase  shift  r~’log£>,,  the  target  pattern 
(4.18)  is  independent  of  the  initial  conditions  \fr(0,  x).  It  depends  only  on  T, 
on  the  largest  eigenvalue  and  on  its  eigenfunction  0,(x).  Moreover, 
a(x)  =  0  except  in  the  area  |  x  —  x0  j<  R.  If  R  is  small,  then  w,  and  ^  are 
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explicitly  given  by 

«,  =  Dk2, 

2 ttD  +  T  j a(x)log  |  x  —  x0 1  dx 

k{  =  2  exp  —  •  y  4 - •  +  •••,  (4.22) 

T  j a(x)  dx 

<M*)  =  2  ttD  /  |r  Ja(x)dxj  +  0(R2logR)  for  |  x  ~x0|<  R, 

*■(*)  =  *o(*il*"*ol)[l  +  0(  R2  log  R )]  for  |  x  — x0|>  R, 

where  y  is  Euler’s  constant  and  K0  is  the  zero-order  modified  Bessel 
function. 

5.  Interacting  Target  Patterns 

We  now  assume  that  a(x)  ^  0  except  in  the  two  well-separated  regions 
|x  -x,  |<  R,  and  |x  ~x2 1<  R2.  We  also  assume  that  j^R  r«(x,  +  s)ds 
>  0  for  both  i  =  1  and  i  -  2.  We  solve  (2.9)  for  bounded  initial  conditions 
ip(0,  x),  and  we  find  that  \(/(T,  x)  is  given  by  (4.9)— (4. 1 2)  exactly  as  before. 
The  only  change  is  that  some  of  the  discrete  eigenfunctions  4>,(x)  are 
centered  around  x  =  x ,  and  satisfy 

%(x)  ~  CJ(6x)(kJri)~x/Ze~k'r'  forr,»l,  (5.1) 

where  /•,  and  0,  are  defined  by  |  x  -x,  j  =  r,(cos  0,,  sin  0, ).  The  other 
discrete  eigenfunctions  are  centered  around  x  =  x2  and  satisfy 

%(x)  ~  Cy(62){k/2y'/2e-k>r>  for  r2  »  1.  (5.2) 

where  |x-x2|=  r2(cos02,  sin  02).  In  (5.1)  and  (5.2),  kj  =  \juj/D  as  before. 

Let  w,  be  the  largest  eigenvalue  whose  eigenfunction  4>,(x)  is  centered 
around  x  —  x,.  Similarly,  let  w2  be  the  largest  eigenvalue  whose  eigenfunc¬ 
tion  4>2(x)  is  centered  around  x  =  x2,  and  assume  that  w,  >  w2.  Then  in  the 
region  where 

e«lr-*1|T-r1|<K  !  e^iT  k2\x-x2\^,  i  r>  i  (5  3) 

the  H(T,  x )  term  in  (4.10)  dominates.  Thus,  A(t,  x)  is  given  by  (4.14)  as 
before.  So,  in  this  region  the  system  evolves  into  a  spatially  constant  bulk 
oscillation. 
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In  the  region  where 

j  ^  eu1T-k1\x-x1\  j  (5.4) 

the  first  term  in  (4.10)  is  much  larger  than  all  the  other  terms.  So  here 

A(r,  x)  =  B([l  +  €<o,/r ]r  +  r_l log fe,4>,(f))  +  0(c),  (5.5) 

analogous  to  (4.18).  Additionally,  when  \x  —  x,  |»  I,  the  argument  of  B 
in  (5.5)  reduces  to  [1  +  cto,/r ]r  —  k]  |  Jc  —  x,  \/T  +  •  •  •  .  Therefore  (5.5) 
represents  a  target  pattern  centered  near  x  =  jc,.  Similarly,  in  the  region 

eo.2T-k2\T-x2\^,  I  g^r-ki ^  (5,6) 

we  have 

A(/,  x)  =  B([l  +  cw2/T]/  +  r_llogZ>2*2(jr))  +  0(c),  (5.7) 

which  represents  a  target  pattern  centered  near  x  =  x2- 

Finally,  when  T  is  large  enough  there  is  an  overlap  region  where 

e^,T-kl\x-Ti\>  j  e»iT-k2 (5  8) 

For  this  region  we  keep  the  first  two  terms  in  (4.10),  obtaining 

A(f,x)  =  B  jr  +  r_llog|  2  /^W'T$/(*)J  j  +  0(c).  (5.9) 

This  equation  describes  how  target  patterns  interact. 

At  any  T  and  Jc  where 

(5.10) 

the  j  —  1  term  is  much  larger  than  the  j  =  2  term  in  (5.9).  Consequently,  at 
this  T  and  x  we  can  neglect  the  latter  term,  and  so  (5.9)  reduces  to  (5.5). 
Thus,  at  this  T  the  target  pattern  centered  at  jc,  extends  to  the  point  x. 
Similarly,  at  any  T  and  x  where 

e«2r-*2|r-i*2|>  ew,r-*,ir-x,i>  ^  (5.11) 

(5.9)  reduces  to  (5.7).  So,  the  point  x  is  part  of  the  target  pattern  centered  at 
x2  at  this  T.  Since  u,  >  w2  has  been  assumed,  at  any  point  x  eventually 

(5.10)  will  hold.  Thus,  as  T  increases  the  target  pattern  centered  near  x  =  x, 
encroaches  upon  the  other  target  pattern.  Eventually  the  second  target 
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pattern  is  engulfed;  it  will  be  unobservable  when 

(5.12) 

If  a(x)=0  except  in  n  >  2  well-separated  areas  i  = 

1,. . .  ,n,  an  analysis  like  the  preceding  one  can  be  done.  Around  each  point 
xt  there  is  a  target  pattern  with  frequency  (1  +  tul/F)P  1  occupying  the 
region  where  »  l.  As  T  increases,  the  domains  of  neighboring 

target  patterns  will  overlap.  Then  the  target  patterns  with  the  larger  w’s  will 
encroach  upon,  and  eventually  engulf,  their  neighbors.  When  T  is  large 
enough,  only  the  target  pattern  with  the  largest  w  will  remain. 


6.  Extensions 

The  analysis  in  Sections  4  and  5  pertains  only  to  two  spatial  dimensions. 
An  identical  analysis  can  be  used  for  one  spatial  dimension.  The  only 
change  is  that  the  (k^y 1/2  factors  in  (4.4),  (4.19),  (5.1),  and  (5.2)  are 
absent  in  one  dimension.  However,  the  analysis  for  three  dimensions  hinges 
on  whether  (4.2)  has  a  discrete  eigenpair  $,(*)  with  «,  >  0.  If  Ta(x)  is 
large  enough  over  a  large  enough  volume,  then  such  a  discrete  eigenpair 
would  exist.  If  such  an  eigenpair  exists,  then  the  analysis  and  results  are 
virtually  identical  to  the  one  and  two  dimensional  cases.  Otherwise  only 
H(T ,  x)  and  the  term  analogous  to  »(7\  x)  would  be  present  in  (4.10),  and 
no  spherical  target  patterns  would  develop. 

The  analysis  can  also  be  extended  readily  to  include  systems  with  delayed 
reaction  terms,  which  often  occur  in  population  models.  For  example,  (1.1) 
can  be  replaced  by 

A,  =  f|a  ,J  -  s,  x))  +  (D„  V2A 

+«g|x,A,jf  tf(s,A(/ -  s,  x))  ds  j.  (6.1) 

If  (6.1)  has  a  stable  periodic  solution  A (t,x)  =  B(/)  =  B(r  H-  P)  at  t  =  0, 
then  we  again  obtain  (2.9)— (2. 11),  and  the  analysis  proceeds  as  before. 


7.  Remarks 

We  used  perturbation  techniques  to  derive  (2.9)— (2.1 1)  from  (1.1),  and 
then  solved  (2.9)  to  obtain  target  patterns.  By  using  perturbation  theory,  we 
have  restricted  attention  to  solutions  which  are  near  the  stable  limit  cycle  B 
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for  all  t  and  x.  There  may  be  other  stable  solutions  of  (1.1)  which  are  far 
from  the  limit  cycle  at  some  x  for  all  t  (10]. 

Also,  the  perturbation  scheme  we  used  tacitly  assumes  that  g(x.  A)  varies 
significantly  only  on  the  0(1)  and  larger  length  scales.  This  restriction  is 
unnecessary.  If  g(x,A)  varies  significantly  on  very  small  length  scales,  a 
similar  perturbation  scheme  can  be  employed.  This  alternative  perturbation 
scheme  yields  solutions  of  (1.1)  given  by  (2.9)— (2.  II),  exactly  as  before. 

Additionally,  the  length  scale  in  (1.1)  has  been  chosen  so  that  the 
diffusion  term  is  the  same  order  as  the  inhomogeneous  term  g(x,  A).  This 
was  necessary.  If  we  used  a  shorter  length  scale,  then  the  diffusion  term 
would  dominate  the  inhomogeneous  term.  To  leading  order,  we  would 
obtain  (2.1 1)  and  (2.9)  with  a(x)  =  0.  As  shown  in  Section  3,  if  a(x)  =  0 
and  if  \p( 0,  x)  is  bounded,  then  $(T,  x)  -*  constant  as  T  becomes  large.  That 
is,  after  an  initial  period  of  time,  \ p  would  become  constant  (to  leading 
order)  on  all  length  scales  smaller  than  the  one  we  have  used.  On  the  other 
hand,  if  a  larger  length  scale  was  used,  then  the  cg(x,A)  term  would 
dominate  the  diffusion  term.  This  results  in  an  inconsistent  perturbation 
scheme. 

The  target  pattern  solutions  obtained  in  Sections  4-6  depend  heavily  on 
the  sign  of  I\  If  T  >  0,  then  the  target  patterns  are  outgoing,  their  frequency 
is  higher  than  the  frequency  of  the  bulk  oscillation,  higher  frequency  target 
patterns  engulf  lower  frequency  ones,  and  higher  frequency  target  patterns 
have  higher  wave  numbers.  If  T  <  0,  the  opposite  of  each  of  these  occurs. 
From  experimental  observations  it  is  clear  that  T  >  0  for  both  the  B-Z 
reaction  and  the  aggregation  of  D.  discoideum. 

In  both  the  B-Z  reaction  and  in  D.  discoideum,  the  periodic  solutions  are 
relaxation-oscillation  type  limit  cycles.  In  this  type  of  cycle,  the  reaction 
evolves  slowly  down  a  quasi-steady  branch,  it  then  jumps  to  a  second 
quasi-steady  branch,  it  evolves  slowly  up  this  branch,  and  then  completes 
the  cycle  by  jumping  back  to  the  first  branch.  As  the  reaction  nears  one  of 
the  jump  points,  a  very  small  perturbation  can  cause  the  reaction  to  jump 
prematurely.  So,  even  very  small  inhomogeneities  g(x,  A)  can  give  rise  to 
significant  frequency  shifts  a(x).  Thus  target  patterns  may  be  more  preva¬ 
lent  for  relaxation-oscillations  than  for  other  types  of  periodic  solutions. 

Finally,  in  (12,  14,  15],  a  phase-diffusion  theory  was  developed.  In  this 
theory,  target  patterns  are  given  by  (2.1 1),  where  \p  evolves  according  to 

\f/T  =  D  V2^  +  a(x)  (7.1) 

instead  of  (2.9).  However,  in  two  dimensions  the  solution  of  (7.1)  is 

«M7\  x)  =  H{T,  x)  +  (4nD)- '/a(y)£,(|  x -y  \2/4DT)  dy,  (7.2) 
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where  H  is  the  solution  of  HT  =  D  V  2H  with  H{ 0,  x)  =  ^(0,  x),  and  £,  is 
the  exponential  integral.  Equation  (7.2)  predicts  that  target  patterns  will  not 
engulf  each  other,  that  the  frequency  of  a  target  pattern  will  essentially  be 
the  frequency  of  the  bulk  oscillation,  that  the  region  occupied  by  a  target 
pattern  grows  like  |  x  —  x0 2/57,  and  that  the  distance  between  succes¬ 
sive  rings  of  a  target  pattern  grows  exponentially  in  |  x  —  x0 1 .  We  conclude 
that  target  patterns  cannot  be  explained  by  phase-diffusion  theory. 


APPENDIX  I 


Here  we  use  Theorems  3  and  5  in  [18]  to  estimate  the  size  of  v(T,  x), 
which  is  defined  in  (4.12).  To  do  this,  we  assume  that  a(x)  has  continu¬ 
ous  second  derivatives  and  we  recall  that  for  some  R,  a(x)  =  0  for  all 
|  Jc -x0|>  R.  Also,  since  H(T,x)  is  the  solution  of  a  heat  equation  with 
smooth  bounded  initial  conditions,  we  note  that  there  are  constants  A/0,  A/,, 
and  M2  such  that 

|//(7,x)|<M0,  |  V//(7\x)|<A/„ 

|  v2//(7\x)|<  A/2,  for  all  |  x  —  x0 1<  R,  alir>0.  (A.l) 
First,  we  define 


o(T,  k)  x)a(x)H(T ,  x)  dx,  (A. 2) 


where  *  denotes  complex  conjugate.  Then  at  any  T,  Bessel’s  inequality 
shows  that 


l©(7\ 


<  T2  lim 

V-*  oo 


/  f‘ 

yA(V)  A) 


- Dk 2 


so(T  —  s,  k)  ds 


dk,  (A. 3) 


where  II  II 2  is  the  £2  norm  and  K(N)  is  the  domain  \/N  <\k\s  N. 
Using  the  Cauchy-Schwartz  and  Bessel  inequalities  now  yields 


||c(7’,  x)|||  <  T2  lim  t[T\  [  o(T  —  s,  k)o*(T  —  s,  k)  dk 

N-oo  J o 

<  r2r JT\\H{s,  x)a(x)||2  ds. 

From  (A.l)  we  now  obtain 

lli>(7\  x)\\2<\T\M0\\a(x)\\2T. 


ds 

(A.4) 

(A.5) 
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Since  D  v  2[H(T ,  x)a(x)]  +  Ta2(x)H(T ,  x)  is  also  in  £2,  a  sequence  of 
steps  like  (A.2)-(A.5)  shows  that 

\\dv2v(t,  x)  +  r<x(x)v(T,  x)ii2<|r|[A/0|r|ii«2n2 

+  Z){M0||v2a||2  +  2M,  II  |  V a  |  II 2  +  M2\\a\\2}]T.  (A.6) 

So  there  is  a  constant  MA  such  that  ||  v  2v(T,  x)||2  <  A/4jT  for  all  7">:  0. 
Thus,  v  and  V  2v  are  both  in  £2,  so  a  Sobelev-type  inequality  shows  that 

|  v(T,  x)  |<  ]/2n  (\\v(T,  x)||2  +  II V  2t>(7\  x)||2)'/2.  (A.7) 

Therefore,  there  is  a  constant  C  such  that 

\v(T,x)\<CT  for  all  x  and  all  T  >  0.  (A.8) 

For  our  purposes,  (A.8)  is  an  adequate  estimate  for  v  in  the  region 
|  x  —  x0 1  <  R.  We  now  obtain  sharper  estimates  for  v  in  the  region  |  x  —  x0 1 
>  R.  By  substituting  (4.10)  into  (4.1),  we  discover  that  v  is  the  solution  of 

m 

vT  -  D  v  2v  =  ra(x)o  +  Ya(x)H(T,  x)  —  2  aj(T)$j(x),  (A.9) 

7=1 


v(0,x)  =  0, 


where 


Oj(T)  =  rf*j(x)a(x)H(T,x)dx^f*}(x)dx  '.  (A.10) 


Therefore,  we  apply  the  Green’s  function 

G(T-s,x  -y)  =[4ttZ)( T  —  5 )] _ '  exp 


4  D(T- s) 


to  (A.9)  and  obtain 


t>(7\ x) |<|y°|  +|/' |  +  2  \Jj\. 


(All) 


(A. 12) 


where 


J°  =  T  fT  f  G{T  -  s,  x  -y)a(y)v(s,  y)dyds , 
■'o  Ja 

Jx  =  r  fT  f  G{T  -  5,  X  -y)a(y)H(s,  y)  dyds, 
Jo  Ja 

Jj  =  fT j  G(T- s,x~y)oJ(s)QJ(y)dyds. 


(A. 13) 


(A. 14) 


Here  Q  is  the  region  |  x  —  x0 1  <  R. 


(A. 15) 
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Equation  (A.8)  shows  that  j  a(x)v(T,  x)  |<  BXT  for  all  T  >  0,  for  some 
constant  Bv  Therefore 

|J0|<B,|r|7'f7'/' G(T  —  s,  x  ~y)  dyds,  (A. 16) 

•'o  Ja 

M°l<fi,|r|r^c-<'iog(i  +  i),  (a.17) 

where  q  =  (\x  —  ,r0 1  —R )2/4DT.  Similarly,  j  a( x)H(  T,  x)  j ^  B2  for  all  T  S 
0,  for  some  constant  Bv  So 

|J>  B2|  r|  ^^-<'1°g(!  +  i).  (A.I8) 

Finally,  there  is  some  constant  Nj  such  that  j  o,(7')$/(x)  |<  for 

all  x  and  all  T  >  0.  Thus 


\Jj\<NjJTJg(T—  s,  x  —  y)e  *'|v  Xnldyds ,  (A. 19) 

N  OC 

| /!<—  f  re~/‘>r~(r~,x~x°l)2/4Drlog{  1  +4DT/(r-\x~xQ\)2}dr. 

2  D  Jq 

(A. 20) 


From  (A.20)  we  find  that 

e-\7-?a\2/A DT  if  kj  |  X  -A0  |<  2 UjT, 

euiT~ki\*~xn\  t  if  kj  |  x  —  x0 1>  2 UjT. 

(A.21) 

When  |  x  —  x0 1  —  R  »  J4DT,  we  use  (A.  12),  (A.17),  (A.  18),  and  (A.21)  to 
estimate  v(T,  x).  Otherwise  we  use  (A.8).  In  particular,  we  note  that  there 
are  positive  constants  Q  and  C,  such  that 

fi{eu'T<bx(x)  >Q[l  +  (A.22) 

and 


1  +  2- 


1  +  2 


1*  ~*ol 

/ DT 

I*  ~*ol 

(DT 


H(T ,  x)  >  C„ 

for  all  a  and  all  7>0.  Thus  |  v(T,  x)  j«  Pxeu,TQx(x)  +  H(T,  x)  when 
T  »  1. 
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